Possible 7r-phase shift at interface of two pnictides with antiphase s-wave pairing 
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We examine the nature of Josephson junction between two identical Fe-pnictides with anti-phase 
s-wave pairing. 7r-phase shift is found if the junction barrier is thick and the two Fe-pnictides are 
oriented in certain directions relative to the interface. Our theory provides a possible explanation 
for the observed half integer flux quantum transitions in a niobium/polycrystal NdFeAsO loop, and 
attributes the 7r-phase shift to intergrain junctions of Fe-pnictides. 

PACS numbers: 



Since the discovery of a new class of unconventional 
superconductors based on Fe compounds, their pairing 
symmetry has been one of the most interesting issues. 
Different from the high Tc cuprates, where the pairing 
symmetry is universally d-wave, there are experimental 
evidences that pairing symmetry in the iron pnictides 
may not be universal. The superconductors with higher 
transition temperature T c are supported by spin singlet s- 
wavei~— , while LaFePO with lower T c seems to have nodal 
in its gap function-. It is interesting to note that the 
iron pnictide has both hole and electron Fermi pockets, 
and the predicted s-wave pairing state has superconduct- 
ing order parameters with opposite signs on the electron 
and hole pockets, often called anti-phase s-wave or s±- 
wave stated—. Among the experiments in support of the 
s± pairing state, the phase sensitive experiment reported 
by Chen et alA provides the most convincing evidence, 
where they observed integer and half integer flux quan- 
tum transitions in a niobium/polycrystal NdFeAsO loop. 
The observed half integer flux quantum demonstrates the 
existence of 7r-junctions in the loop of niobium and poly- 
crystal pnictides, hence a direct evidence for the opposite 
signs of the superconducting order parameters in different 
Fermi pockets. In passing, we recall that phase sensitive 
experiments provided a direct evidence for the d x 2_ y 2 
pairing in superconducting cuprate s 12 ' 13 . 

Because of the polycrystal nature in the sample of the 
NdFeAsO, a phase shift of ir in the composite loop in 
tunneling could occur at the Nb-Nd-1111 interface, or at 
the junction between two Nd-1111 grains 4 . Theoretically, 
there have been several studies to examine the possible tt- 
phase shift involving an interface between a conventional 
superconductor and a Fe-pnictide with s± pairing under 
certain conditions^—, which may help understand the 
possible 7r-junction at the Nb-Nd-1111 interface. The 
possibility of a 7r-phase shift in the interface of two Fe- 
pnictide intcrgrains of the same doping has not been care- 
fully examined, although intuitively one may consider it 
unlikely. In this paper, we study a Josephson junction 
between two iron pnictides. When the junction barrier 
is thick and the two Fe-pnictide grains are oriented in 
certain directions relative to the interface, the junction 
could give a 7r-phase shift. 

The 7r-junction is a Josephson junction with negative 



critical current, and the critical current for junction be- 
tween conventional superconductors is defined as 
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where Tk q is the tunneling matrix, £j(k) = 
^€i(k) 2 + Aj(k) 2 is the quasiparticle energy of the su- 
perconductor i = 1, 2 respectively, £j(k) is the single elec- 
tron energy measured relative to the chemical potential, 
and Ajm) are superconducting gap functions, which we 
shall assume to be real here. The co-efficient of the inte- 
gral will always be taken to be positive throughout this 
paper. Eq. ([T]) can be easily generalized to the iron pnic- 
tide superconductors with multi-bands, and the critical 
current are given by, 
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where a and j3 are index of bands, and all other notations 
are the same as in the single band expression except for 
the additional band indices. As pointed out by Sigrist 
and Riceii, though the critical current is not gauge in- 
variant, the parity of the number of the 7r-junctions in a 
loop is gauge invariant. So in the following, we choose 
a convenient gauge where the gap functions of the hole 
pockets are positive and the gap functions of the electron 
pockets are negative. 

In the usual case, the junction between two identical 
superconductors is a 0-junction with I c to be positive. 
Let us consider a point junction of two identical pnictide 
superconductors. In this case, = to, independent of 
crystal momentum and the band. The critical current is 
found to be 
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According to Ambegaokar and Baratoff, the above for- 
mula can be further written as 
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where N a is the density of states of Fermi pocket a, 
A m = min(|A Q |, |A^|) and A M = max(|A Q |, |A^|) are 
the smaller and larger gap on the two pockets respec- 
tively, and K is the elliptical integral. In the spe- 
cial case, all the gap functions have the same ampli- 
tude, the elliptical function K(0) = n/2, so that I c oc 
(22 sgn(A a )N a ) 2 > 0. This positive definiteness ap- 
pears to remain valid when the gap amplitudes are dif- 
ferent. In particular we consider below iron based su- 
perconductor in the 2-dimension limit. We shall work in 
extended Brillouin zone for convenience and set the Fe-Fe 
distance to be length unit. There are two electron Fermi 
pockets with one around (ir, 0) and one around (0, 7r), and 
two hole Fermi pockets around (0,0) point. We consider 
the case, where the superconducting gaps on one of the 
hole pockets and on the electron pockets are about the 
same while the gap on another hole pocket is smaller, as 
reported in ARPES for BaFe2 As^. By using the proper- 
ties of the elliptical function, Vl - k 2 K(k) < K(0) = § , 
the critical current I c is found to be always positive. This 
illustrates that the different signs of the s— wave gap func- 
tions is necessary condition, but not sufficient condition 
for a 7r-junction. Below we shall examine a thick barrier 
junction between two pnictide superconductors of certain 
orientation and show that such a junction may give a tt- 
phase shift. 

We consider two half-infinite iron pnictide samples sep- 
arated by a vacuum barrier with height U and width d 
as shown in fig.[T] The interface is along yz plane so that 
the momentum components along y and z directions are 
conserved in the tunneling process. For a thick barrier, 
the usual assumption that the tunneling matrix element 
T is independent of momentum or crystal momentum is 
no longer valid. This can be illustrated by examining the 
free electron tunneling process in quasi one-dimension, 
where |T| 2 is the transmission coefficient of the scatter- 
ing problem. For a potential barrier normal to x direction 
with height U and length d, the transmission coefficient 
reads 

ijni2 _ ^K 2 k x q x 

~ K 2 (k x + q x ) 2 + (k 2 + k 2 x ){n 2 + ql) sinh 2 ( K d) ' 

(5) 

where k x and q x are the ^-direction wavevector of the 
incoming and outgoing plane waves respectively, k is 
the imaginary wavevector inside the barrier with k = 

j-./2m(U — E) + h 2 k 2 , m is the mass of the electron, 



and fc|| is the wavevector parallel to the barrier interface, 
which is conserved in the scattering process. In the thick 
barrier limit, i.e. nd 3> 1, one have 

\T\ 2 1& * 2 K<lx 2 Kd (6) 

If we assume that p = j- ^2m(U — E) ^> ku, it can be 
further simplified as 
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so the transmission coefficient decays exponentially with 
the increment of planar wave- vector k\\ . The above for- 
mula can be extended to the electron tunneling process 
in a lattice with the following modifications as pointed 
out by Mazing. All the wave-vectors in the expression 
of the above equation are replaced by the correspond- 
ing group velocities except k\\ in the exponential factors 
which tracks the oscillation of the wavefunction paral- 
lel to the interface direction. The second modification 
is to replace the plane wavefunction in the free electron 
case by a Bloch wave ip n k(r) — e lk ' r u;„k( r )- The peri- 
odical function w n k can be further Fourier's transformed 
as w„k(r) = J2k ^nk,Ke' K r , where K is reciprocal lat- 
tice vector. If u;„k(r) is localized, one can approximate 
Fnk.K. ~ constant for not very large K. The wavefunc- 
tion can finally be written as 



^k(r)=F, ik 5>* (k+K) - r . 
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According to eqn. l[Tj). the tunneling is mainly from the 
component with K\\ = 0, so the tunneling matrix reads 

l T kq I 2 « ^ak,/3q ^^^(k+K^Xq+QJn 
K,Q 



(h 2 n 2 + m 2 v 2 akx )(h 2 n 2 + m 2 v 2 0clx ) 



(9) 



where K and Q are the reciprocal lattice vec- 
tor of FeAs sample 1 and 2 respectively, k — 
\^j2m(U - E) + H 2 (k + K) 2 , the delta function tracks 

the planar momentum conservation, A is a factor related 
to the detail information of the electron wavefunction 
and can usually be approximated as constant, and v a u. x 
and vpq X are the group velocity along x direction of the 
electron in band a(f3) with lattice wavevector k(q) re- 
spectively. 

The exponential factor in above equation can be 

_i d < 2k+K 'ir K n d , , n 
rewritten as e * e p , where the first term is 

the contribution of the Brillouin zones with K\\ = and 

the second term is the additional factor for the Brillouin 

Zones with finite Ku. In the following, we ignore the z- 

axis dispersion , which is very weak in 1111 compound, 

and only consider the physics of FeAs plane. For the 

configuration shown in fig. [T] we have that K y = 2nn/a 

for FeAs sample 1 and Q y = \f2m/Kja for FeAs sample 2, 

where a is the lattice constant. For a barrier with d ~ 10a 

lOir/a, the largest factor for 



and p = \yj2m(U - E) 
Brillouin zone with finite K\\ is e~ 27r « 0.002. So only 
the contribution from the Brillouin Zone with Ku — is 
important. And because p 2 ^> k 2 , eqn. (0 becomes 
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Because of the fast drop of the coherence factor A 



2E 



in the critical current when the state is away from the 



3 




FIG. 1: Schematic diagram of a thick Josephson junction of 
two iron pnictides separated by a vacuum. The Fe lattice 
orientation for the two materials are depicted in the upper 
panel. Lower panel is the corresponding Fermi surface in the 
Brillouin zone for the two material. 

Fermi surface, the planar momentum conservation paral- 
lel to the interface, and the exponential factor in tunnel- 
ing matrix in eqn. (|10[) , only the electron pockets around 
X and the hole pockets of iron pnictide sample 1 and the 
hole pockets of sample 2 are important in the configura- 
tion shown in Fig. [TJ Because k|| ~ 0, the exponential 
factor in eqn. (1101) becomes a constant. For simplicity, 
we assume that the gap A is momentum independent 
which is consistent with the ARPES result^. By substi- 
tuting eqn. (fT0| into eqn. ([2]) and noticing that velocity 
v x oc -§£-■> the critical current reads 

Er r e M( k v) r e M( k v) A^A' 3 
/ dk y I dei / de 2 wr — - w-, 

(11) 

where e^(fc y ) and e^(fc a ) are the minimum and maxi- 
mum energy of the electron in band a with given k y and 

any k x . If e™ (/3) (fcj,) - E f < -A a()3 ) and e^\k y ) - 
Ef 3> A. a (p-\, the integral can be approximated with 

sgn(A Q A^)A m X( v /l - A^/A^), where A m and A M 
are the smaller one and larger one of amplitude of the 
two gaps A Q and A^. If the amplitude of the supercon- 
ducting gap on those bands are very close to each other, 
the elliptical integral K can be approximated as a con- 
stant, and the final result is 

I c oc ^2 sgn(A ct A /3 )A m min(Ai a , A2/3), (12) 

a/3 

where A; 7 is the width along the direction parallel to the 
interface of the Fermi pocket 7 of FeAs sample i. 

Then we consider the real materials. Most iron-based 
superconductors in 1111 family are electron doped in- 
cluding the material used in the experiment by Chen et 
ali^ The theoretical calculations and experiments show 
that the parent compound have two hole pockets around 




FIG. 2: The type of Josephson junction. The red solid line 
is X hh /X eh = A^/Aj^ 1 which separate the two regions. See 
text for details. 

r point and one electron pockets around X point. Af- 
ter doping of the electrons, one would expect the elec- 
tron pocket expands and the hole pockets shrink, when 
the electron concentration is large enough, the inner hole 
pocket will be too small to contribute to the critical cur- 
rent or even vanished. In this case, eqn. fjl 2|) becomes 

I c oc A^X hh — A^X eh , (13) 

where A^f = min(|A lQ |, |A 2/3 |), = mm(X la ,X 2p ). 
So the Josephson junction is a 7r-junction when 
X hh /X eh < A^/A^ and vise versa, as shown in fig. [2] 
The ARPES experiment has shown that the gaps on elec- 
tron pockets and the outer hole pocket are same 2 -, i.e. 
Ajj' 1 = A^. So the sign of the critical current is de- 
termined by the width of the Fermi pockets along the 
direction parallel to the interface. The DFT result of 
LaFeAsO has suggest that the hole pocket has a little 
anisotropy and the Fermi wavevector along (110) direc- 
tion is a little larger than the one along (100) direction 20 , 
i.e. X'i < X' 2 l and X hh — Aj. On the other hand, if the 
electron concentration is large enough, one should have 
that the width of the electron pocket is larger than the 
one of the hole pocket, i.e. X\ > X 2 and X eh — X 2 > X hh , 
which indicates that the junction is a 7r-junction. An- 
other possibility is that the electron concentration of 
sample 1 is a littler larger than the one of sample 2. In 
this case, the size of hole pocket of sample 1 should be 
smaller than the one of sample 2, i.e. X\ < X 2 . Again, if 
the electron concentration is large enough, one should 
have Af > X 2 , and the junction should also be a ir- 
junction. 

For the material in 122 family, both electron and hole 
can be doped with various ways. For the electron doped 
122 material, the ARPES experiment suggests that there 
is only one hole pocket around T point, and the widths of 
hole pocket and electron pocket are comparable 19 which 
should be similar with the e-doped 1111 case. So one 
should observe similar phenomena in the e-doped 122 
material with large enough electron doping. But the sit- 
uation is different in the hole doped case where both of 
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the two hole pockets will contribute to the critical cur- 
rent, and the size of hole pockets should be larger than 
the one of electron pocket. By taking account that the 
gap on outer hole pocket is same as the one on electron 
pocket, one can find that the critical current is positive 
according to eqn. (fT2|) . and one can not observe half- flux 
jump in those hole doped 122 material. 

In summary, we have examined the Josephson junc- 
tions between two Fe-pnictides with the same electron 
concentration. We found that it is impossible for two 



Fe-pnictides to form a 7r-junction in the case that tun- 
neling matrix element is momentum independent. But a 
7r-junction can be formed in the thick barrier case where 
the pockets with large parallel momentum are irrelevant. 
Such 7r-junctions may be the origin of the 7r-flux jump in 
C.-T. Chen et al.'s experiment. 
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